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Abstract 

We consider the thermal equilibrium distribution at inverse temperature f3, or 
canonical ensemble, of the wave function ^ of a quantum system. Since I? spaces 

00 ' 

■ contain more nondifferentiable than differentiable functions, and since the thermal 

equilibrium distribution is very spread-out, one might expect that ^ has proba- 
bility zero to be differentiable. However, we show that for relevant Hamiltonians 
I/" - ) \ the contrary is the case: with probability one, ^ is infinitely often differentiable 

and even analytic. We also show that with probability one, ^ lies in the domain 
of the Hamiltonian. 
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1 Introduction 

We address the question whether the wave function \1/ of a typical system from the 
canonical ensemble of thermodynamics with inverse temperature f3 is differentiable. As 
pointed out in [7j, the thermal equilibrium distribution of the wave function, correspond- 
ing to the canonical ensemble, is the "Gaussian adjusted projected measure" GAP(p), 
a probability measure on the unit sphere in Hilbert space whose definition we recall in 
Section 2, for p = p@, the density matrix of the canonical ensemble, given by (see, e.g., 

181 EH) 

pp = z(jhrf~ PH with ^^) = tre_/3H - (!) 
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Thus, we take \l/ to be a random unit vector with distribution GAP(pp). The surpris- 
ing result is that in many relevant cases ^ has probability one to be infinitely often 
differentiable and even analytic, i.e., GAP{p p ){C°°) = GAP{pp){C w ) = 1. 

We explore four kinds of arguments concerning the smoothness of each requiring 
somewhat different assumptions on the Hamiltonian H and leading to somewhat differ- 
ent conclusions. Some of the arguments do not depend on the special measure GAP(pp) 
but show that, for suitable Hamiltonians H, every distribution whose density matrix is 
Pp will be concentrated on the smooth (resp. analytic) functions; other arguments use 
the way the measure GAP(p) is constructed from a Gaussian measure. The measure 
GAP(p) is discussed in detail in [7j. It has density matrix p and is stationary if p is. 

The first argument aims at showing that the Fourier coefficients of ^ go to zero so 
fast that they are still square-summable after multiplication by any power of the wave 
number k. This can be easily applied to cases in which the eigenf unctions of H are 
plane waves, such as for the free Schrodinger equation in a box. The second argument 
is based on the assumption that the eigenf unctions of H are smooth, and the theorem 
asserting that, for a series of functions, summation and differentiation commute if the 
series of the derivatives converges uniformly. We formulate a condition on H that entails 
this kind of convergence almost surely (a.s.) for the expansion of \& in eigenfunctions 
of H. The third argument, which supposes that \l/ is a function on an interval I C R, 
aims at showing that the increments are not too large, + Aq) — < Aq, which 

suggests differentiability; however, the rigorous version of this argument provides only a 
very weak result. The fourth argument, the simplest and most elegant one, is of a more 
abstract nature: it concerns not the question ^ G C°° but instead the related question 
\l/ e domain(if^) for f e N; indeed, we obtain without further assumptions on H that 
a.s. \1/ G C°°(H) := P)fc=i domain (if ) for almost all (3 for which thermal equilibrium 
exists at all. We also provide a variant of this argument concerning the space of analytic 
vectors of H. 

The paper is organized as follows. In Section 2, we recall from j7] the definition of 
the measure GAP(pp) representing the canonical ensemble. In Section 3, we study as 
an example of H the Laplacian on the circle. We conclude smoothness and analyticity 
of \1/ from an analysis of the decay behavior of the Fourier coefficients of In Section 4, 
we apply the same argument to the (relativistic or nonrelativistic) ideal gas in a box. In 
Section 5, we take into account the symmetrization of the wave function for describing 
bosons or fermions. In Section 6, we give the second kind of argument, providing a 
general criterion on the Hamiltonian that is sufficient for concluding that \1/ is a.s. 
smooth. The criterion concerns bounds on the derivatives of the eigenfunctions of H . 
In Section 7, we discuss the third argument, which concerns the direct estimation of the 
difference quotients of ^ . In Section 8, we describe the fourth argument, which allows 
to conclude that \l/ lies in the domain of H and all its powers. In Section 9, we conclude, 
as an application of our results, that \1/ a.s. possesses a Bohmian velocity field. 



2 



2 The Canonical Ensemble 



In this section we give the definition of the measure GAP(p) on the unit sphere y(Jtf) 
of Hilbert space as introduced in [7j. 

The measure GAP(p) is defined for every density matrix p (positive operator with 
tr p = 1) on Jjf. We obtain the thermal equilibrium measure GAP(pp) by using the 
canonical density matrix (|TJ) for a self-adjoint operator H (the Hamiltonian) and a 
number (3 > (the inverse temperature) such that 

Z(J3,H) = tre- pH < oo. (2) 

The measure GAP(p) is defined as the distribution of the random vector 

^ GAP = ^ GA /||^ GA ||, (3) 

where \I/ Gj4 is a random vector with distribution GA(p) (the "Gaussian adjusted mea- 
sure") defined by 

GA{p){&i,) = U\\ 2 G{p)m i (4) 

where G(p) is the Gaussian measure on M 3 with covariance matrix p. 

More explicitly, for a random vector ^ G to be G(p)-distributed means that for any 
0i,02 G the components Zi = (0i | ^ G ) and Z 2 = (0 2 |* G ) of # G are complex 
Gaussian random variables with mean zero and covariance 

EZxZ* = E(0i|^ G )(^ G |0 2 ) = (0i|p|0 2 ), (5) 

where E denotes expectation. In particular, if {|v? n )} is an ortho normal basis of Jtf 
consisting of eigenvectors of p with eigenvalues p n , then the coefficients (ip n \^/ G ) of \I/ G 
are independent complex Gaussian random variables with mean zero and variances 

E\( Vn \q G )\ 2 = Pn . (6) 

If p is of the form Q then the ip n are also eigenvectors of H. 

Although we are interested only in those p of the form ((TJ for physically relevant 
Hamiltonians, we will sometimes, when this makes the mathematics clearer and more 
elegant, formulate facts or conditions in terms of an arbitrary density matrix p. 

For any probability measure p on y(Jif), its density matrix is given by 

P,= j MdVOMM, (7) 

or Pn = E /t |-0)(?/>|, where E^ denotes the expectation with respect to p. In particular, 
E^l^l^) = ((plp^Lp) for every fixed ip G M 3 . (If a probability measure p on is 
not concentrated on y(J$?), the notion of density matrix of p does not make sense any 
more; however, (J7J), or E M |^>)(^|, is still the covariance matrix of p.) As mathematically 
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expressed by E^^P^) = tr(p /i P) for every projection P, p M provides the distribution 
of outcomes of any quantum experiment on a system with /i-distributed random wave 
function. The density matrix of GAP(p) is indeed p. This and other fundamental 
properties of the measure GAP(p) are discussed in [7]. 

The following simple fact will sometimes be useful as it reduces the task of showing 
smoothness of ^ = ty GAP to showing smoothness of the Gaussian random vector \I/ G 
with distribution G(p). For any subspace W of Jif, we have that 

if G(p)(W) = 1 then GAP(p)(W) = 1 . (8) 

To see this, note that GA(p) has the same null sets as G(p) (as it is absolutely contin- 
uous with respect to G(p) with a density that vanishes only at one point), so that, if 
G(p)(W) = 1, = G(p)(Jf \ W) = GA(p)(JT \ W) and thus GA(p)(W) = 1; but, by 
definition ©, if ^ GA G W then also ^ GAP G W. 



3 A Case Study: The Laplacian on the Circle 



H = ~ — ^ (9) 



In this section we consider a single particle moving on a circle S 1 with the free Schrodinger 
Hamiltonian, 

2m dq 1 

where m denotes the mass of the particle, q the angular coordinate on the circle, and 
wave functions are written as periodic functions of q. The result we derive is that 
relative to any measure p on S^(Jif) with density matrix pp with (3 > (or, in fact, any 
measure p on Jtif with covariance matrix pp), almost every wave function ip is smooth, 
ip G C°°(R); we then go on to show that //-almost every wave function is analytic, 
i) G C^iR). 

We begin with considering, instead of differentiability, a closely related property: ex- 
istence (in L 2 ) of the distributional derivative. In other words, we consider the property 
of a wave function tp that \k\ ip{k) is still square integrable where ip is the Fourier trans- 
form of if}. Since the functions we are considering are 27r-periodic in q, the appropriate 
property is that the Fourier coefficients c^, defined by 

#?) = ]Tc fe e ifc o, (10) 

fcez 

are still square-summable after multiplication with \k\. Let W denote the £-th Sobolev 
space, i.e., the subspace of L 2 ([0,27r]) containing those functions whose Fourier coeffi- 
cients Ck satisfy 

\k\ e Ck is square-summable, i.e., \k\ 2e \cu\ 2 < oo. (11) 
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We ask whether ip G W for a random wave function ip with distribution \i. Since the 
eigenf unctions of H are the plane waves, 

^,(g) = -Le^, neZ, (12) 

V Z7T 

the energy coefficients of a wave function are just the Fourier coefficients. The eigenval- 
ues are 

h 2 

E n = T^n 2 , neZ. (13) 
Thus, our question about ip amounts to asking for which I G N we have 

^^|(^|^)| 2 <oo. (14) 

This indeed holds /i-a.s. for all t G N; to see this, note that 1 

E^ri*\(<p n \i(>)\ 2 = J> 2 %|<^)| 2 = (15a) 

1 \- f^-^L # —1— V n * e -W#l»nW < ^ (15b) 



because for any constant 7 > and for sufficiently large |n|, 

(2£ + 2) log Inl < 7|ra| 2 and thus \n\ 2t e' 1 ^ < (16) 

\n\ 2 

If the expectation (|15|) of a [0, oo]-valued random variable is finite, the variable is a.s. 
finite. Thus, 

p(Z>*M^>r = 00) = or fi(W £ ) = 1 (17) 

for all t G N. 

We now make the connection between and C £ , i.e., with classical differentiabil- 
ity: by the Sobolev lemma p. 52], every function in W e is equal Lebesgue-almost- 
everywhere to a function in C l ~ l . Hence, every function in f]^i W e is equal Lebesgue- 
almost-everywhere to a function in C°° . In particular, /i-a.s. there is a G C°° such 
that ^(9) = 0(9) Lebesgue-almost-everywhere. 

We now turn to analyticity. By a similar argument as given in (|15|) and one 
can see that 

/i(^e 2Q H|^ n |^)| 2 = 00 ) =0 (18) 

for every a > 0, so that 

e Q ' fc ! Cfc is a.s. square-summable. (19) 



1 Numbers on top of equality signs indicate which equation is being applied. 
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Regarding the variable q in (fTUj) as complex, we observe that the right hand side of (fTUj) 
converges, as a consequence of (fTU|). uniformly in every strip —a + e < Imq < a — e 
with < £ < a. (To see this, we use that square-summable sequences are bounded, 
e a ' fc ' |cfc| < C, so that for q in the strip, \cke lkq \ = \ck\e~ klmq < \ck\ e^ a ~^ < Ce'^ £ , 
which is summable over k £ Z.) Since the uniform limit of analytic functions on an 
open set in the complex plane is analytic (by virtue of the Cauchy integral formula), ip 
is analytic in the strip —a < Img < a; since a was arbitrary, ip is entire (i.e., analytic 
on the whole complex plane). More precisely, p-a.s. there is an entire function cf) such 
that ip(q) = <p(q) Lebesgue-almost-everywhere in KL 



4 The Ideal Gas in a Box 

In a similar way, we can treat any Hamiltonian whose eigenfunctions are plane waves. 
A particularly relevant case is that of the ideal gas: N noninteracting particles in a 
c?-dimensional box [0,7r] d with Hamiltonian 

N ft 2 

tf = -£-A, (20) 

1=1 

with Dirichlet boundary conditions, where Aj is the Laplacian acting on the coordinates 
of the i-th particle. Our conclusion will again be that p(C°°) = p{C w ) = 1 for every 
measure \i on whose covariance matrix is p@, and in particular for p = GAP(pp). For 
the moment, we ignore the symmetrization postulate; we will treat bosons and fermions 
in Section 5. 

The Hamiltonian H on = L 2 ([0, ir] Nd ) has eigenfunctions 3, p. 78] 

/ 2 \ Nd/2 JL. JL 

<Pn(q)={-) l|||smKa9i, a ) (21) 

i=l a=l 

where n = (ni ( i, . . . , n N ^ d ) G N Nd and q = . . . , qN,d) £ [0, n) Nd , and eigenvalues 

n d 2 2 
E n = YY— n 2 ia = — INI 2 - (22) 

i=l a=l 

The right hand side of ()21j) extends in an obvious way to a function on M. Nd that is 
27r-periodic in every variable, which we also call (p n ; using the coefficients ((p n \4>) of the 
energy expansion we have a natural extension of any ip e ffl to a function <fi on "§L Nd 
that is 27r-periodic in every variable, = ^ n (y>n|^) Vn- Its Fourier coefficients are 

N d 

c k = (-^) 7vd (nn si g n (M)(^i fcl , 1 i,..,i^ d i|^) (23) 

i=l a=l 
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where k = (k^x, . . . , fcjv )£ j) G "L Nd and we set sign(O) = 0. 

We begin with the existence in L 2 of the £-fold distributional derivative. We assert 
that /i-a.s. for all I e N 

J] ||A;|| 2£ | Cfe | 2 <oo. (24) 

kez Nd 

This follows from 

X) ll^iri^| 2 = 2 M £ ||n|| 2 %|<^>| 2 = (25a) 

m O^d imn O^d 

® Z V llnll^e-^ ~ V ||„||2^ p -(/3?i 2 /2m)||n|| 2 < McU 

- z m m 2^ ll n ll e - z{(5 H) 2^ W n \\ e ^ ^ 5b J 

< E ( 1 + "v)* ■ ■ ■ (! + *- m ,2m)M = (25c) 

^' ' n€N Nd 

= ( 2 I> + v?' e -^^) m t =o, (25d) 

where we used 



I ll / j fc ; t^j ■ ■ \ - 7 — fc^Lfc/ 



Inequality (|2*4"j) means that <fi lies in the Sobolev space W , and by the Sobolev lemma 
[3 p. 52] also in C m for all m < i — Nd/2. Since I was arbitrary, is /x-a.s. smooth, 
and thus so is ^, its restriction to [0, 71"]^. 

The same argument can be applied to the relativistic case, in which the Hamiltonian 
is the free Dirac operator 



N 



H = - ^ i icha i ■ v * + ^c 2 A) (27) 

i=l 

with c the speed of light, m the mass, and a, and $ the Dirac alpha and beta matrices 
acting on the i-th spin index of the wave function. Again, one obtains that /i(C°°) = 1 
for all fi with covariance matrix pp. 

We turn to analyticity and to this end assert that /i-a.s. for all a > 

£ e 2Q ll fe " |c fc | 2 < oo. (28) 

kez Nd 

This follows from the fact that, by the same reasoning as in (|23|) . 

£ e2a " fc " l^l 2 = zWm E e 2a Nl e -(^/ 2m )NI 2 <00 
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because for any constant 7 > and for all but finitely many n e N Nd , 2a||n|| — 7|M| 2 < 
— (7/2) ||n|| 2 , while e - *- 7 / 2 -"'™" is summable over N Nd by By the same argument as 
in the last paragraph of Section 3, one can conclude from (J28|) that <fi is analytic in the 
cylinder {q G C Nd : ||Img|| < a}. Since a was arbitrary, <p is entire. Thus, /i-a.s. there 
is an entire function such that ip(q) = <p(q) Lebesgue-almost-everywhere in [0, 7r] . 
(For the Dirac equation, since the energy eigenvalues grow like c/i||fc||, ip a.s. possesses 
an analytic continuation to the cylinder ||Img|| < f3cH/2.) 

5 Bosons and Fermions 

In the previous section, we ignored the symmetrization of the wave function for systems 
of bosons or fermions. If one takes the symmetrization into account, one reaches the same 
conclusion: smoothness is almost sure. But instead of going through the calculation 
of the previous section again, we provide a simple argument why GAP(pp) must be 
concentrated on C°° for indistinguishable particles (with symmetrized wave functions) 
if it is concentrated on C°° for distinguishable particles (with unsymmetrized wave 
functions). 

The symmetric (respectively anti-symmetric) state vectors form a subspace of Jif = 
L 2 ([0,ii} Nd ); let P denote the projection to that subspace; the subspace can be written 
PJf. Since the Hamiltonian (j2"Uj) is invariant under permutations, we have that HP = 
PH = PHP. Thus, the canonical density matrix for indistinguishable particles is 

P(i(PH, PJ^) = c Ppp{H, J^)P (29) 

where we have made explicit the dependence of pp on the given Hamiltonian and Hilbert 
space, and c = Z{(3, H)/Z(f3, PHP). 

Now observe that for a Gaussian measure with covariance p, we have that G(PpP) = 
G(p) o P _1 , where P~ l is understood as mapping subsets of PJtf? to their pre-images in 

, in other words ^ G{ - p P p ^> = P^/ G ^ in distribution. 

In our case, P is the symmetrization (respectively anti-symmetrization) operator, 
which maps smooth functions to smooth functions. Since is a.s. smooth (by the result 
of the previous section), so is P^ G ; with © we conclude that GAP(pp(PH, PJ^))(C°°) = 
1. The same argument works with analyticity. 

6 A General Sufficient Condition for Smoothness 

We now present a second kind of argument, different from the one used in the previous 
sections; it applies to the measure GAP(p) but not to all measures with density matrix 
p. The argument provides us with a condition, see (J3Uj) below, on any given density 
matrix p (and thus, for p = p@, on the Hamiltonian) ensuring that GAP(p)(C°°) = 1. 



Theorem 1 Let Q be an open subset of M. d . Suppose that the density matrix p on 
= L 2 (Q,C m ) has C°° eigenfunctions ip n (q) with \\<p n \\ = 1 and eigenvalues p n , such 
that for all n and all £ = 0, 1, 2, 3, . . the £-th derivative of <p n is bounded, 



|vV 



„i loo = sup |V <p n (q)\ < oo, 



where by absolute values we mean 



|vV(?)l s 



E E 



(<?) 



l j^\dq h ---dq 



If for all £ = 0,1,2, 3,..., 



^ ||VVn||oo \fp~n < OO 



(30) 



then GAP(p)(C 00 {Q,C m )) = 1- 



Proof. To begin with, for a complex Gaussian random variable Z with EZ = and 
ElZl 2 = a 2 one can determine that 



E\Z\ 



dx dy — — exp 



7T(J Z 



x 2 + y 2 



7T 



-a. 



(31) 



Setting Z = (ip n \^ G ), we obtain that 

E E|| V Vn(^n|^ G ) =Ej2\(<Pn\V G }\ ||VVr. 



^E|(^ n |^ G )| ||VVn||oo ^ E"V"V^II V Vn||oo "< OO, 



and therefore 



< oo = 1. 



Prob(^||vVn (<Pn\V G ) 

n 

Since this is true of every £, we have that in the expansion 



(32) 



(having C°° partial sums), a.s. the £-th derivatives of the partial sums converge uni- 
formly; and in particular, (|32j) itself converges uniformly. It is a standard theorem 
(see, e.g., jH p. 118]) that if a sequence f n of functions converges pointwise and the 
derivatives V/ n uniformly, then the limit function / is different iable and has derivative 
V/ = lim V/ n . Therefore, a.s. \I/ G G C°°(Q), and the derivatives are 



(33) 



By (JEJ), a.s. \1> G C°°(Q), which completes the proof. □ 

By applying this proof to local coordinates, we can generalize the result to Rieman- 
nian manifolds and vector bundles as follows. Let Q be a Riemannian C°° manifold, 
E a C°° complex vector bundle over Q with positive- definite C°° Hermitian inner prod- 
ucts on the fiber spaces, and let V be the covariant derivative operator corresponding 
to a C°° connection on E. Let M 3 = L 2 (E) be the Hilbert space of square-integrable 
(with respect to the Riemannian volume) measurable cross-sections of E, and C°°(E) 
the space of smooth cross-sections. Suppose that the density matrix p on Jif has C°° 
eigen-cross-sections (p n (q) with \\(p n \\ = 1 and eigenvalues p n , such that for all n and all 
£ = 0, 1, 2, 3, . . the £-th covariant derivative of <p n is bounded, 

||VVn||oo = SUp |VVn(g)| < OO, (34) 
<?eQ 

where the absolute values are taken with respect to the Riemannian inner product on tan- 
gent spaces and the Hermitian inner product on fiber spaces. If for all I = 0, 1, 2, 3, . . 

^ ||VVn||oo y/P^ < OO (35) 

n 

then GAP(p)(C 00 {E)) = 1. 

Another easy generalization of Theorem Q concerns analyticity: Let Qc be an open 
subset of C d and Q := {(qi, ■ ■ ■ ,qd) G Qc '■ ?i> • • • > Qd G ^} Q ^- d - Suppose that the 
density matrix p on Jif = L 2 (Q,C m ) has eigenvalues p n with normalized eigenvectors 
ip n G C w (Qc); with C u (Qc) the space of L 2 functions on Q that possess analytic con- 
tinuations to Qc- If f or every compact set K C Q c , 

J^IM^IL VP™ < °°> ( 36 ) 

n 

(writing also (f n for the analytic continuation and ^p u \k f or Us restriction to K), then 
GAP(p)(C«(Qc)) = l. 

Proof. By flgfl and (|3T|). 

e 2JW k IL K^»l* G )l = z2\\^\k\\ 00 ^Vp^< oo, 

n n 

and thus a.s. ^ n ||<£n|A:|| [(v 9 "!^ )! < oo. As a consequence, the expansion 

* G (q) = J2vn(q)(<Pn\* G ) (37) 

n 

converges not only for q G Q but also for q G Qc, in fact uniformly on every compact 
set K C Q c . Since uniform limits of analytic functions are analytic, \1/ G and thus also 
ty GAP are a.s. analytic. □ 
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In order to demonstrate that the conditions (J3*U|) . and are not unreasonably 
strong, we show that they are satisfied for the Laplacian on the circle. Here the eigen- 
functions are given by (|12j) . and their derivatives, respectively analytic continuations to 
disks K = {q G Qc = C : \q\ < a}, have bounds 

||VVn||oo = -i, llynklL = ^=e a|n| . 



/2tt 00 V2tt 

In this case, ()30|) respectively f)35|) is satisfied since 



Eii*vgi-«- iw - = 4=E 



2ir 



U|/ e -09ftV4m> 2 <00 



by (JEJ), and similarly ()36)1 since 



'2-7T 

7 Estimating Difference Quotients 

In this section, we follow another line of reasoning for studying regularity properties of 
\P, based on a standard theorem on the regularity of sample paths of a Gaussian process. 
However, the result is much weaker than what we obtained in the previous section. 

Assume for simplicity that the configuration space is an open interval ICR, and that 
Jif = L 2 (I, C). The idea in this section is to consider the increments \1/ G (g + Aq) — ^ G (q) 
of the Gaussian process \I/ G and to argue that for reasonable Hamiltonians they are of 
the order of magnitude of Aq > 0, 

\V G (q + Aq)-V G (q)\ < Aq, (38) 

which suggests that difference quotients converge to differential quotients as Aq — > 0, 
i.e., that \& be differentiable. However, what can rigorously be concluded from a state- 
ment about the variance of the increment analogous to (}3*H|) is less than differentiability, 
namely Holder continuity with exponent 1 — e. 

We now describe the argument in detail. We pretend that \1/ G is everywhere defined 
in I (although, strictly speaking, vectors in Hilbert space are equivalence classes of 
functions coinciding Lebesgue-almost-everywhere) in such a way that it is a Gaussian 
process in the sense that, for any choice of qi,..-,q n £ I, the joint distribution of 
\1/ G (gi), . . . , \E fG (g n ) in C n is Gaussian. It follows that for any Aq > 0, the increment 
ty°(q + Aq) — \I/ G (g) is a Gaussian variable, and we can compute its variance, 

E|^ G (g + Aq) - ^ G (q)\ 2 = p(q + Aq,q + Aq) - 2Re p(q, q + Aq) + p(q, q), (39) 

where p(q,q') = (q\p\q') are the "matrix elements" in the position representation of the 
density matrix p = pp. Assuming that 

p(q, q') is a smooth function, (40) 
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which would appear to be a reasonable assumption on the Hamiltonian, we can employ 
a Taylor expansion of p to the second order around (q, q) and obtain from that 



E|^ G (g + Aq) - ^ G (q) 



Gi 



d 2 p 



dq dq' 



Aq 2 + 0(Aq 3 ). 



(41) 



q'=q 



It is a standard result jHl Thm. 8 of Chap. Ill] that for a Gaussian process \1/ G with the 
following bound on the variances of the increments: 



E\^ G (q + Aq) - ^ G (q)\ 2 < K Aq p , 

where K > and p > are constants, the realization a.s. satisfies 

|^ G (g + Aq) - ^ G (q)\ < K' Aq p/2 | log Aq\ 1+S 

for arbitrary 5 > and a suitable constant K' = K'(5) > 0. Inserting (}4*T]) into 
obtain from (JI3J), in case d 2 p/dqdq' ^ 0, that a.s. 

\^ G (q + Aq) - ^ G (q) \ < K" Aq 1 ' 6 

for arbitrary e > 0, i.e., Holder continuity of degree 1 — e. 

In order to obtain a stronger estimate than (|44j) . one might hope that 



(42) 



(43) 



we 



(44) 



d 2 p 



dq dq' 



for all q E I. 



(45) 



q'=q 



This would allow us to replace the exponent in (|44|) by 3/2 — e, which would give us 
in particular local Lipschitz continuity, so that ^ would be differentiable Lebesgue- 
almost-e very where. But any exponent greater than 1 is too good to be true. Indeed, 
as we shall show presently, (J43j) holds only for one particular density matrix po, the 
projection to the 1-dimensional subspace of constant functions. For p = p Q , ^f(q) is 
constant with modulus determined by normalization and random phase. 

To see p = p Q , write p in terms of its eigenfunctions f n (q) and eigenvalues p n , 
P(Q,Q') = Y, n PnPn(q) <(<?')> and observe that 



d 2 p 



dq dq' 



q'=q 



q'=q 



where ip' n denotes the derivative of ip n . The only way how this quantity can vanish for 
all q is that all (p n have identically vanishing derivative and thus are constant, which 
implies p = p . 
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8 Concentration on the Domain of H 



In this section we utilize a fourth kind of argument, different from those of the previous 
sections. It is our most elegant argument, particularly simple and direct, as it deals only 
with the eigenvalues and eigenvectors, but also more abstract. This argument applies 
not only to GAP(pp) but to every probability measure p, on S^{J^) with density matrix 

The question we address is whether p, (domain (H)) = 1. The answer is yes. Even 
more, we show that, for any self-adjoint H and almost all /3 for which Z(/3, H) < oo, any 
p, with density matrix pp is supported by the domain of H for every I G N. (Whether 
■0 G domain(if) implies differentiability depends of course on H.) After that, we show 
further that GAP(pp) for sufficiently large (3 is concentrated on the subspace of analytic 
vectors pQ of H. 

Theorem 2 Let p be a density matrix on the Hilbert space Jff, p, a probability measure 
on with density matrix p, and f : [0, oo) — > K a measurable function. If 

tr(pf(p) 2 ) <oo (46) 

then p, (domain (/(/?))) = 1, where the domain of f(p) can be defined, in terms of an 
orthonormal basis {\ip n )} of eigenvectors of p with eigenvalues p n , by 

domain(/(p)) = \jp G Jrf? : ^|/(p„)(^ n |^) | 2 < oo}. (47) 

n 

Proof. From (|46j) it follows that 

I / (P«) <V»« I ^> 1 2 =J2f(Pn) 2 ^ J ( i Pn\^)\ 2 = ^f(Pn) 2 Pn < OO. 

n n n 

Therefore, /x-a.s. Zln|/(P™)(v 5 n|^) | 2 < °°> or A t (domain(/(p))) = 1. □ 

Corollary 1 Let H be a self-adjoint operator on the Hilbert space ffl" ' . Suppose Z((3q, H) < 
oo for some (3q > 0, which implies Z(f3, H) < oo for every (3 > /3o- Then for every /? > /3 
and every probability measure p, on with density matrix pp, pl(C°°(H)) = 1, where 

C°°{H) = ^ domain^). 

Proof. For p given by (0), we have that H — — -| logp + i? id for some constant E . 
Define f(x) = (— ^logx + E ) e for x > and f(x) = for x = 0. Since f(p) = H e , 
Theorem El yields the claim if we can confirm the condition ([46)1. which we do now. 

Since tr exp(— (3qH) < oo, there is a basis {|<p«) : n G N} of eigenvectors of H with 
eigenvalues E n . Furthermore, only finitely many of the eigenvalues lie below zero, so 
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that the set M := {n G N : E n > 0} contains all except finitely many numbers. Observe 
that for every (3 > /3 , 



oo > tre-*> H = e ~ P ° En > Yl e ~' 30En > XX 

neN n&N n£Af 



■PE n 



and thus Z((3, H) < oo. For e > with e < (3 — /3 , we find, for any £ G N, 

oo > ^ e -(/3-e)^ = J- e"^ e^ > ^ e"^ X ^ = X |[ X E * e"^. 

n&Af neJV raeAf fc=0 ' fc=0 ' neAT 

In particular, 

^^e-^<oo, 

n&A/" 

which implies tr (ppH 2 ^ < oo. □ 

If P > 2/3o (with Z([3 ,H) < oo), we obtain the stronger result for the measure 
GAP(pp) that it is concentrated on the subspace C U (H) of analytic vectors of H, i.e., 
those vectors ip G C°°(H) with 

t^<~ (48) 

for some £>0|lj. It is sufficient for if) £ C w (if) that 

JV 1 ^ 1 |<n|0}| < oo (49) 



because then 

oo 



oo _^ oo 



^=0 £=0 n 



< 



oo 



n £=0 ' n 

For < £ < /3/2 - /3 , flUJ) is a.s. true of t/> = \1/ G , and thus also of -0 = \I /G ' AP , because, 
assuming without loss of generality that all E n > 0, we have by (|31|) that 



E V e £E " I I = V e £E " V^_ e -^/ 2 < ^E__ y 
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9 Existence of Bohmian Velocities 



As a final remark, we mention an application of smoothness of the wave function: 
differentiability is needed in Bohmian mechanics [2], a theory ascribing trajectories 
to the particles of nonrelativistic quantum mechanics. This is because the Bohmian 
law of motion, which for N particles with masses mi,...,mjv at the configuration 
Q(t) = (Q 1 (t),...,Q N (t)) reads 

-7T = vJ(Q) = — Im— — Q , 
dt rrii ip*ip 

involves the derivative of the wave function. Suppose the wave function ip is chosen 
at random according to the canonical distribution GAP(pp) with inverse temperature 
(3. Then any condition on the Hamiltonian entailing that ip is a.s. smooth also implies 
that the Bohmian velocity vector field on configuration space Q = M 3iV , whose i-th 
component is vf , is a.s. well defined everywhere outside the nodes of ip. 

The analogous conclusion holds, as we shall explain presently, for the numerous 
further variants of Bohmian mechanics that have been considered (such as Bohmian 
mechanics on curved spaces, on the configuration space of a variable number of particles 
jH], for wave functions that are cross-sections of a complex vector bundle, and variants 
suitable for the Dirac equation or for photons). The laws of motion of these variants, 

are defined by giving the appropriate expression for the velocity vector field on the 
manifold Q, and these definitions of xr can be summarized by the formula jHj 

4)., /(!) = B /!?») V/ecns). (50) 

Here, / : Q — > R is an arbitrary smooth function with compact support playing the 
role of a coordinate function, and numerator and denominator involve inner products 
in the value space of ip (which may be a fiber space of a vector bundle of which ip is 
a cross-section). For if) G domain(if) and / G C£°(Q), the right hand side of (joTIj) will 
be well defined since multiplication by / maps the domain of H to itself, since H is 
the sum of a differential operator (of up to second order) and a multiplication operator. 
Since the /'s from C^°(Q) suffice for determining (up to changes on a null set), one 
obtains indeed a vector field v^, defined on Q \ {q : ip(q) = 0}, for every ip from the 
domain of H. 

Hence, every wave function ip from the domain of H is sufficiently regular to define 
a Bohm-type velocity field. By Corollary ^ the random wave function \1/ with the 
thermal equilibrium distribution GAP^pp) possesses a velocity field v * with probability 
one, provided that there is (3q < (3 with Z(f3 , H) < oo. 
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